Neutron beta decay and the current 
determination of V u d* 



A. Garcia, J.L. Garcia-Luna 1 and G. Lopez Castro 2 

Departamento Fisica, Centro de Investigation y de Estudios Avanzados del IPN, 
Apdo. Postal 14-740, 07000 Mexico, D.F., Mexico 



Abstract 

Measurements of neutron beta decay observables required to determine \V U( j\ are 
reaching the 0.1% accuracy. In this paper we review the calculation of the decay 
rate of this process, discuss its relevant uncertainties, and obtain an expression that 
is precise at the 10 -4 level. Our analysis clearly shows the necessity of more precise 
measurements of A, the ratio of axial/vector couplings. The current situation in 
neutron beta decay is that one cannot yet quote a single consistent value for \V u d\ 
from it. We also discuss the region of parameter space in the IK^-A plane where 
new physics effects should lie, if they contribute to neutron beta decay. 



1. Introduction. Neutron beta decay can provide in the future the most 
precise determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix el- 
ement \V u d\ [1]. The experimental inputs required for this purpose are the 
neutron lifetime r n and the ratio of axial-vector to vector couplings, A = 
(?i(0)//i(0). An update of the precise computation of the neutron decay rate 
is important for (at least) the following reasons: (a) the measurement of r n 
and A have reached the 0.1% accuracy and further improvements are expected 
in the future (see for example [2]) and (b) independent and precise determi- 
nation of V u d is compelling to test the unitarity condition J2i=d,s,b I Km 1 2 = 1 
of the CKM matrix. This requires from theory the computation of the decay 
rate at the 10~ 4 level in accuracy, and will allow a determination of \V u d\ with 
a precision 0.1% which becomes competitive with the determination coming 
from superallowed Fermi transitions (SFT) among J p = + nuclei [3], namely, 

\V ud \ = 0.9740 ±0.0010 . (1) 

The purpose of the present paper is to review the computation of the neutron 
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beta decay rate and discuss its relevant uncertainties at the 10~ 4 level. At this 
level it becomes necessary to discuss the relevance of the induced hadronic 
form factors and even radiative corrections up to a 2 contributions [4,5]. 

There are two reasons why jV^I can be obtained from SFT with very small 
error bars. One is of statistical nature, because Eq.(l) arises from the weighted 
average of 9 measured transitions [6]. The other one is an advantage of SFT 
over neutron beta decay; they involve at tree- level only the vector current. The 
present uncertainty in \ V u d\ from SFT has two dominant (theoretical) sources: 
the nuclear-dependent (mismatch) isospin breaking corrections and the piece 
of photonic corrections induced by the axial- vector current [7,8]. The advan- 
tage of neutron beta decay is that it is free of a nuclear environment, but 
its (theoretical) disadvantage is that it involves at tree-level the axial-vector 
current, too. Away from a nuclear environment, the isospin-breaking correc- 
tions are more reliable to compute. Such corrections clearly affect both leading 
form factors. However, the isospin breaking corrections enter into /i(0) only 
at second order in the breaking by virtue of the Ademollo-Gatto theorem [9]. 
This fact makes them very small and they contribute only at the 10~ 5 level 
[10]. On the other hand, these corrections are more important in <7i(0) and, 
also, they are more difficult to compute reliably. This complication is avoided 
using the experimental value of <?i(0), through the measurement of the ratio 
A. The main source for such measurements is the electron-neutron-spin asym- 
metry coefficient. It is because the experimental uncertainty in A can be much 
improved over its theoretical uncertainty that the neutron beta decay deter- 
mination of \V u d\ may become its best determination. 

For sakeness of clarity we will present the different ingredients for the deter- 
mination of the decay rate in the following order. In section 2 we provide the 
expression for the decay rate of the neutron at tree-level by including small 
effects of non-leading form factors (other than /i(0) and gi(0)). In section 3 
we address the issue of radiative corrections. We will discuss the uncertainties 
associated with factorization of different orders in a as commonly presented 
in the literature. In section 4 we summarize the current information of ex- 
perimental data on r n and A including the most recent results. In particular, 
we will discuss the discrepancies in different measurements of A, and we will 
provide a comparison of the different determinations of \V u d\ as obtained from 
different sources. Section 5 contains our conclusions. 

2. Decay rate at the tree- level. The tree-level amplitude (mean- 
ing that radiative corrections are ignored) for the decay n — > pe~u e can be 
parametrized in terms of six Lorentz-invariant form factors. The most general 
form of this amplitude is given by (we follow the conventions of Ref. [5]): 
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where G is the effective weak coupling of the relevant four-fermion hamilto- 
nian at the tree-level and Z M is the leptonic V-A current. ^ denote dimen- 
sionless q 2 = (p — p') 2 dependent form factors. Since the momentum transfer 
is very small (qmax/ m n ~ 1-9 x 10~ 6 ), we can already anticipate that the 
contributions of the form factors /2,/3,<72 and g 3 to the decay rate are very 
suppressed. By the same reason, the q 2 dependence of form factors can be 
safely neglected. 

By including recoil effects (in the phase space factor, too), the rate of neu- 
tron beta decay can be written as [5]: 
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where m e is the mass of the electron and /j, gi are form factors evaluated at 
zero momentum transfer. The coefficients A/ 4 , X 9i denote the corrections due 
to q 2 dependence of /$ and We do not display the contributions of f 2 , g 2 , 
Xf 2 , \ 92 and \ g3 since their coefficients are at the 10~ 8 level, i.e., well below 
the current experimental precision. In Eq. (3) we improved the precision of the 
coefficients of the leading form factors by one more decimal place, with respect 
to Ref. [5]. 

In order to get a simpler expression that incorporates the terms in Eq. (3) 
relevant at order 10 -4 let us discuss the size of the different contributions. As 
is well known, in the limit of exact isospin symmetry the conserved current 
hypothesis (CVC) allows one to obtain /f yc = 1 and f? vc = 1.8529 [5]. Also, 
exact G-parity, that is expected to hold in the isospin symmetry limit, predicts 
the absence of second class form factors: fs = and g2 = 0. The contribution 
of the axial form factor g 3 to the decay rate is expected to be below the 1CT 4 
level, even if, according to the PCAC hypothesis, g 3 = (2m 2 l /m 2 )gi 114. 

When isospin symmetry gets broken (due to the mass difference of u and d 
quarks), the above predictions change. For example, /i(0) will deviate slightly 
from unity while (?2 and f 3 become different from zero. An estimate of these 
effects indicates that 1 - f x w 5 x 1(T 5 [10], g 2 ^ Q x 10~ 3 and f 3 ~ 10" 3 [11]. 
Therefore, the expression for the decay rate at the lowest order that is precise 
at the 10~ 4 level gets the simple form: 

i? = ^|K,| 2 (1.6335)(l + 3A 2 ). (4) 



3. Radiative corrections. An expression for the decay rate that is 
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precise at this level requires the inclusion of radiative corrections up to order 
a 2 . The formal perturbative expansion in powers of a for the decay rate can 
be expressed as follows: 

R = R° + R(a) + R(a 2 ) + ■■■ (5) 



where R° is the rate at tree-level discussed in the previous section and R(a n ) 
denote its corrections of order a n . 

At any given order in a, the radiative corrections can be split into two 
terms: R(a n ) = Rui^a 11 ) + i?MZ)(a n ), the model-independent (MI) and the 
model-dependent (MD) pieces, following their independence or not of the de- 
tails of weak and strong interactions. 

The MI radiative corrections, also called inner corrections, arise from long- 
distance QED (virtual and bremsstrahlung) corrections up to an energy scale 
of 0(1 GeV). The MD radiative corrections, also called outer corrections, in- 
volve essentially the short- distance electroweak corrections. 

In the following sections we will discuss in detail the leading orders of ra- 
diative corrections that have been computed in the literature. 
3.1 Corrections of order a. The corrections of order a have been computed 
long ago [12,13]. The model-independent piece of 0(a) was given as [12] 

R M i(a) = ^i|VU| 2 (l + 3A 2 )^ J (f a + ^9{E)) dS , (6) 



while the model-dependent contribution as [13] 
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In the above equations E denotes the energy of the electron, I — \l\ its three- 
momentum and dS = El(E m — E) 2 dE . The function f a denotes the coulombic 
correction of order a, 
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and the expression for g(E) can be found in Ref. [12]. The origin of the different 
terms in Eq. (7) can be found elsewhere [13,7]. Numerically, the MI radiative 
corrections at this order are 

— b I f a dS = 0.0546 , (9) 
ml J 
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As we can observe, the dominant piece of MI radiative corrections of order 
a comes from the coulombic term. We will leave the numerical evaluation of 
the MD corrections to the following subsection. 

Before we close this subsection let us emphasize that the remaining effects 
of radiative corrections that have a universal nature have been absorbed in the 
redefinition of the effective weak coupling constant with the result Go —> Gp, 
where Gp = 1.16639(2) x 1CT 5 GeV -2 can be extracted from muon decay [3]. 
3.2 Corrections of order a 2 . In this subsection we consider the MI radiative 
corrections of 0(a 2 ) [14]. The MD corrections of this and higher orders will 
be discussed in subsection 3.3. In a similar way to the previous subsection, we 
can write the MI corrections as [4]: 



where again, we have separated explicitly the coulombic correction f a 2 from 
the remaining corrections 5 a 2 at this order. 
Numerically, we have: 



As in the case of order a, we observe the dominance of the coulombic piece in 
the MI radiative corrections of order a 2 . 

3.3 Corrections of order higher than a 2 . The leading pieces of MI ra- 
diative corrections of order a 3 have been computed in Ref. [14]. Following 
the pattern already observed, we can expect that the MI corrections will be 
dominated also by the coulombic term. Indeed, coulombic corrections can be 
included to all orders (via the Fermi function F(Z,E) [15]) by solving the 
non-relativistic Schrodinger equation for the electron in the Coulomb field of 
the proton. However, it it clear that for the purposes of our analysis, we have 
to keep only the MI corrections up to 0(a 2 ). 

Now, we will focus on the MD radiative corrections. Following the pro- 
cedure discussed in [13,7], let us rewrite the MD corrections of 0(a) in the 
form 
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The term ln(mz/m p ) is by far the dominant piece of this correction. As in 
the case of coulombic corrections, one may wonder about the relative size of 
higher order MD corrections having logarithmic terms. Ref. [7] has pointed 
out that the effects of higher order can be approximated by resummation of 
the leading-logarithmic corrections of 0(a n \v?{m z /m p )), n — 1, 2, • • • via the 
renormalization group equation (RGE). The net effect of this resummation is 
to replace Eq. (14) by [7]: 
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where R' MD denotes the MD radiative corrections that take into account dom- 
inant logarithmic corrections at higher orders. 

According to ref. [7], S(m P) m z ) can be obtained by solving an specific RGE 
subject to the boundary condition condition S(m p ,m p ) = 1. The numerical 
result for this correction is: 



S(m p ,m z ) = 1.0220 . 



(16) 



To obtain this numerical value, which differs slightly from Ref. [7], we have 
used the fact that the top quark is heavier than the Z boson and should not 
be included in the thresholds when solving the RGE. 

3.4 Final remarks on radiative corrections. If we collect the results pre- 
sented in the previous sections, we can write the radiative corrections to the 
neutron decay rate as 
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In some terms of the above equations we have replaced a _1 (0) = 137.036 — > 
a^irrip) = 133.93. According to Ref. [7] m p is the relevant scale for the elec- 
tromagnetic coupling in the renormalization scheme adopted in these calcula- 
tions. 

Let us stress that Eq. (17) has the correct additive structure in terms of the 
expansion parameter a. This observation is very important in view of several 
factorization approximations currently found in the literature [4,8] which can- 
not be justified when the decay rate is computed at order 10~ 4 . For example, 
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it is customary to take the coulombic corrections to all orders through the 
Fermi function as a global factor in the integrand of Eq. (17) [4]. This would 
lead to the following expression for the decay rate 



a(m p ) 
1 + -^9(E) + 5 a 2 
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Clearly, this factorization introduces spurious terms of order a 2 , through the 
product of the term of order a in the Fermi function with other contributions 
of 0(a) in the remaining corrections. A similar argument can be applied to the 
factorization of outer and inner radiative corrections (see ref. [7] and Eq. (7) 
in Ref. [8]). Therefore, in the computation of the decay rate at the 10~ 4 level 
(or to a higher accuracy) radiative corrections must be taken in an additive 
form as shown in Eq.(5) in order to avoid spurious corrections of order a 2 (and 
higher orders). 

After these remarks, we proceed to obtain a numerical expression for the 
decay rate involving radiative corrections. The numerical values of some of 
these corrections have been provided in Eqs. (9,10, 12, 13 and 16). Currently, 
the largest uncertainty in radiative corrections comes from the terms which 
originate in the axial-induced photonic corrections, 
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The non-asymptotic piece C has been re-evaluated in ref. [8] with the re- 
sult C = 0.881 ± 0.030. The mass parameter m A is taken as a low energy 
cutoff applied to the evaluation of the axial part of the asymptotic box 
correction. Marciano and Sirlin [7] have suggested that, given the ambigu- 
ities in the choice of this cutoff, it is more appropriate to take the range 
400 MeV < rriA < 1600 MeV for this axial mass. This wide range, is at 
present the major source of uncertainty in radiative corrections to neutron 
beta decay, since the term given in Eq. (19) becomes 

A = 0.0022 ± 0.0008 . 



If we insert the different values of radiative corrections discussed above, we 
obtain an expression of the decay rate that is accurate at the 10~ 4 level: 

R = 0.1897|Kd| 2 (l + 3A 2 )(1 + 0.0739 ± 0.0008) x 10~ 3 s" 1 , (20) 



where the quoted error comes from the uncertainty in the axial cutoff param- 
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eter [7,8]. 

Notice that if we had included radiative corrections with the Fermi function 
F(Z, E) and MD corrections in a factorized form, we would have obtained, 

R = 0.1897|Kd| 2 (l + 3A 2 )(1 + 0.0758 ± 0.0008) x 10~ 3 s' 1 , (21) 



to be compared to Eq. (20). 

As it was previously mentioned, the factorization of some pieces of radia- 
tive corrections affects the expression for the decay rate at the order 10~ 4 . 
Therefore, in the following we will use our Eq. (20) to quote a determination 

of|Kd|. 

4. Determination of \V u d\. Comparison with other sources. 

According to Eq. (20), two experimental inputs are required to determine \ V u d\ 
from neutron beta decay: the neutron lifetime r n = R~ l and A. We first sum- 
marize here the current situation regarding their measurements and afterwards 
we discuss [V^l- 

4.1 Neutron lifetime. The world average reported by the Particle Data 
Group in the 1998 edition is [3] 

rr° 98 = ^ = (886.7 ± 1.9) s. (22) 



This value is obtained [3] from a weighted average of ten independent mea- 
surements and the quoted error has been inflated by a scale factor S=1.2. A 
more recent and more precise measurement has been reported in the literature 
[6,20,22]: r n = (885.4 ± 1.0) s, which is consistent with Eq. (22). Therefore, 
from its weighted average with Eq. (22), we obtain 



The set of measurements of r n is consistent as reflected by a low value 
of the scale factor (S=1.07). As it was mentioned in the introduction, the 
measurement of r n has reached the impressive accuracy of 0.1 %. 
4.2 Ratio A = <7i(0)//i(0). The consistency of measurements of A is not clear 
at present. So far, measurements of the electron spin asymmetry have provided 
the most precise values reported for A. By combining the four most precise 
measurements (typically those with an accuracy higher than 0.4%) one gets 
the weighted average [3] 

x pdg98 = L2670 ± Q.0035 , (24) 



where the error has been increased by an scale factor of S = 1.9. Therefore, 
the choice of this set of measurements would indicate a precision of 0.3 % in 
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A. Other less precise measurements may be excluded because they have neg- 
ligible statistical weig ht at the 10" 4 level. 

In order to study the impact of different measurements of A on the deter- 
mination of \V u d\, we prefer to quote the value of this parameter as a function 
of different values of A. Following Ref. [3], we will choose the four most precise 
measurements of A: Aj, where i = L,B,Y,R denote, respectively, the results 
of references [18,19,16,20]. We will exclude from our analysis the result of Ref. 
[17], because it has been superseded by the new result in Ref. [20]. The explicit 
values are given in Table 1. The updated average becomes 

x pdg99 = L26g7 ± Q.0016 , (25) 



where the error has not been increased by any scale factor. 

Several comments about Eq. (25) are in order. First, the two most precise 
values of A used to obtain Eq. (25), namely, Ay and A_r differ each one only 
by two standard deviations from Eq.(25) and by 4.7 standard deviations from 
each other. This reflects an inconsistency between both measurements. Second, 
the same inconsistency can be observed in the ideogram shown in p. 622 of 
ref. [3] where the set of 4 measurements used to quote Eq.(24) are grouped 
into two separated peaks of Gaussian distributions. This signals the lack of 
meaning of the single weighted average given in (24) and points out towards 
the necessity of new and independent measurements of A to elucidate the 
nature of the systematic errors. 

4.3 Determination of \V u d\ from neutron beta decay. In Table 1, we show 
the values of \V u d\ obtained for individual measurements of A and some of its 
partial averages. Also displayed are the corresponding values of the associated 
scale factor S = Xi/( n ~ 1) when more than one value of A is used; the 
bigger the scale factor is, the larger the inconsistency within the set of chosen 
data. The inconsistency among the measurements of A can be appreciated in 
several ways (see Table 1): (a) the central values of \V u d\ corresponding to Ay 
and \ R differ by almost 4<r, (b) if we use the two sets of data that appear 
clustered into two separated Gaussian distributions (see section 4.2), namely, 
Xlyb and \r, we find that the corresponding central values of \V u d\ differ by 
more than 4a. 

We observe that in the latter case (b) the corresponding values of \V ud \ 
violate the unitarity of the CKM matrix. Indeed, if we use \V US \ = 0.2196 ± 
0.0023 and |T4b| = 0.0033 ± 0.0008 [3] we obtain (for comparison we include 
the unitarity check for SFT) 



VudV + WusV + W, 



ub 



0.9928 ±0.0031, for \ R , 
0.9969 ± 0.0022, forSFT, 
1.0066 ± 0.0035, for \ LYB 



(26) 
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It is very interesting to observe that when we combine the set of four data, 
namely, if we use Xlbyr, we obtain a value of \V u d\ that is (accidentally!) con- 
sistent with unitarity \V ud \ 2 + \V US \ 2 + \V ub \ 2 = 0.9986 ± 0.0025 . 

Therefore, since the whole set of most precise measurements of A is in- 
consistent, we can only quote two different values of \V u d\, namely, those cor- 
responding to the two clearly separated Gaussian distributions shown in the 
ideogram of p. 622 in Ref. [3]. Explicitly (for comparison we include \ V u d\ from 
the unitarity of CKM matrix) 



IK 



ud\ 



0.9719 ± 0.0019 ± 0.0004 from X R , 

0.9790 ± 0.0016 ± 0.0004 from \ LYB - 

0.9756 ± 0.0005 from unitary of CKM matrix. 



(27) 



The first quoted error arises mainly from the uncertainty in the measure- 
ment of A, while the second one comes from the uncertainty in the calculation 
of MD radiative corrections. This result indicates that besides the interest of 
new measurements of A to solve the current discrepancy, they are also required 
to provide the most precise determination of |T4<i|- 

One may wonder where future measurements of A should lie if no new 
physics effects are relevant in neutron beta decay. In Fig. 1. we have plotted 
the one standard deviation elliptical contours of \V u d\ and A obtained from 
Eq. (20) and the two above values of A. For comparison we also display the 
two horizontal bands corresponding to \V u d\ from unitarity and SFT at one 
standard deviation, too. 

The small ellipse in Fig. 1 within the CKM band is the one standard de- 
viation region allowed by the unitarity \ V u d\ and Eq. (23) and it is the region 
where no new physics can yet appear in neutron beta decay. The current sit- 
uation is that the two values of the \V u d\ from neutron beta decay and its 
value from SFT do not overlap with that region. Thus, Fig. 1. provides a vi- 
sualization of Eqs. (26) and (27). At any rate, neither of these discrepancies 
with unitarity has enough statistical significance to be currently accepted as 
an indication of the existence of new physics. But, the smallness of this region 
illustrates the potential of neutron beta decay to either detect new physics or 
to provide strong lower bounds on its existence. 

5. Conclusions. The determination of \V u d\ from neutron beta decay has 
several advantages. First, isospin breaking effects are well under control since 
they affect the decay rate at the level below the 10~ 4 (contrary to SFT where 
they enter at a few tenths of percent level [8]). Second, the piece of MD ra- 
diative corrections C (see Eq. (19)) which is associated to the non-asymptotic 
axial-induced photonic corrections only receive contributions from one-body 
terms (contrary to SFT where the photon and the W gauge boson appear- 
ing in the box diagram can be coupled to different nucleons). Third, possible 
residual effects of Z-dependent radiative corrections (see Ref. [21]) are absent 
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in neutron beta decay. In addition, the non-leading form factors are irrelevant 
at this level of precision. The complications of determining gi(0) do remain. 
Fortunately, this form factor can be taken from experiment, because its ex- 
perimental value can be made appreciably more precise than its theoretical 
one. Thus, the main source of theoretical uncertainties in the neutron beta de- 
cay rate are the radiative corrections. We have seen that they are well under 
control up to a few parts in 10~ 4 . Prescriptions used to incorporate differ- 
ent contributions of these corrections become relevant. By this we mean that 
the (commonly used) factorization approximations of some pieces of radiative 
corrections introduce spurious contributions at the 10~ 4 level. The current 
prescription is the additive one. 

However, new measurements of A are required, first, to solve the discrep- 
ancy among existing data and, second, to improve the determination of \V u d\- 
Actually (see Eq. (27)), an improvement of a reduction of a factor 2 in the 
error bars of A would provide a determination of | V u d | that is competitive with 
the one coming from SFT. 

The main conclusion in the above analysis is that the theoretical decay 
rate of neutron beta decay is well under control up to a few parts at the 1CT 4 
level. Given the present discrepancy of experimental measurements of A and 
our limitation to combine those data in a statistically significant way, we can 
only quote the two favored values of \V u d\ as given in Eq. (27). These values 
correspond to the two subsets of consistent data according to the ideogram 
shown in p. 622 of the Particle Data Group [3]. New physics effects can only 
be detected in neutron beta decay outside the small ellipse within the CKM 
matrix band in Fig. 1. 

From Eq.(27) it becomes clear that, once the error bars on A are appre- 
ciably reduced, it will be very important to reduce the theoretical error bars. 
Therefore it would be very useful to invest some effort in better controlling 
the model dependent piece of the radiative corrections in neutron beta decay. 
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A 


scale factor 




A L =1.2660(40) 


_ 


0.9766 ± 0.0026 


Ay=1.2594(38) 


— 


0.9809 ± 0.0025 


A B =1.2620(50) 





0.9792 ± 0.0032 


A fl =1.2735(21) 


_ 


0.9719 ±0.0015 


A LB =1.2644(31) 


0.6 


0.9777 ±0.0021 


A Ly =1.2625(28) 


1.2 


0.9789 ± 0.0019 


X YB = 1.2604(30) 


0.4 


0.9802 ± 0.0020 


A Lfl =1.2719(19) 


1.7 


0.9729 ± 0.0014 


Ayfl=1.2702(18) 


3.2 


0.9740 ± 0.0013 


A Bfl =1.2718(19) 


2.1 


0.9730 ± 0.0014 


\ LYB = 1.2624(24) 


1.6 


0.9790 ± 0.0017 


A LB fl=1.2707(17) 


1.8 


0.9737 ±0.0012 


Ay B fl=1.2692(17) 


2.5 


0.9746 ± 0.0012 


Ay Li? =1.2695(17) 


2.4 


0.9744 ± 0.0012 


A L y B fl=1.2687(16) 


2.1 


0.9749 ± 0.0012 



Table 1 

Values of \V U( j\ obtained from different values of A and r = (885.9 ± 0.9) s, via Eq. 
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0.98 - 



0.978 



0.976 - 



. 974 



. 972 - 



0.97 



0.968 









^^^/////////////^^ Matrix ////////// 
















A(LYB) 


, . . , A ( R ), , . . . . " 



1.26 



1.2625 



1.265 



1.2675 

A 



1.27 



1.2725 



1.275 



Fig. 1. \V u d\ vs. A. The two large ellipses corresponds to the la contours obtained 
via Eq. (20). We also display the two horizontal bands corresponding to \V u d\ from 
unitarity of the CKM matrix and SFT at la, too. The small ellipse within the CKM 
matrix band illustrates the potential of neutron beta decay to detect new physics 
or to put strong lower limits to its existence. 
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